Identifying the multiferroic phase of doped CuFe02 using inelastic neutron scattering 
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We report inelastic neutron scattering measurements that provide a distinct dynamical "finger- 
print" for the multiferroic ground state of 3.5% Ga-doped CuFe02- The complex ground state is 
stabilized by the displacement of the oxygen atoms, which are also responsible for the multiferroic 
coupling predicted by Arima. By comparing the observed and calculated spectrum of spin excita- 
tions, we conclude that the magnetic ground state is a distorted screw-type spin configuration that 
requires a mechanism for magnetoelectric coupling different from the generally accepted spin-current 
model. 



Multiferroic materials allow the electric polarization 
to be controlled by switching the direction of magnetic 
ordering and consequently offer prospects for many new 
technological applications [JH] . Because multiferroic be- 
havior has been found in materials that exhibit com- 
plex (non-collinear and incommensurate) magnetic or- 
der, determining the spin arrangement of the ground 
states is essential to understand its role in the multifer- 
roic coupling jiMlfll . In the spin-current model or inverse 
Dzyaloshinskii-Moriya mechanism of multiferroelectric- 
ity j the electric polarization is induced by inversion 
symmetry breaking associated with a complex spin spiral 
structure. This model predicts that the electric polariza- 
tion is perpendicular to both the chiral axis Si x Sj and 
the magnetic ordering wavevector Q. However, in ma- 
terials like CuFeC>2, the electric polarization is parallel 
[ll-13| to Q, indicating that the spin-current model can- 
not explain the multiferroic coupling. The multiferroic 
behavior of CuFe02 can instead be explained 0] by spin- 
orbit coupling in the presence of lattice distortions that 
lower the crystal symmetry. Therefore, it is important 
to understand how lattice distortions affect the magnetic 
structure of doped CuFe02- 

Nakajima et al. 13[ proposed that the magnetic struc- 
ture of multiferroic CuFcC>2 is a "proper" spiral config- 
uration composed of two alternating simple spiral struc- 
tures. Fishman and Okamoto on the other hand, 
reported that the ground state for a frustrated triangular 
lattice has a complex non-collinear (CNC) spin configu- 
ration characterized by the presence of higher spin har- 
monics. Elastic neutron scattering measurements alone 
are not sufficient to distinguish between these two com- 
plex magnetic states because they produce similar sets of 
observable satellite reflections. Consequently, novel ways 
must be devised to determine the true magnetic ground 
state. This Letter reports inelastic neutron scattering 
(INS) measurements that provide a distinct dynamical 
"fingerprint" for the multiferroic ground state of 3.5% 
Ga-doped CuFe02- From our analysis, we conclude that 
the magnetic ground state is a distorted screw-type spin 
configuration that is incompatible with the generally ac- 




FIG. 1: (a) The three-dimensional crystal structure of 
CuFeC>2. (b) The two-dimensional hexagonal lattice struc- 
ture of the Fe 3+ layers and the intralayer interactions. The 
orange and blue colored atoms illustrate the displacements 
of oxygen atoms lying above and below the hexagonal plane. 
The magnetic structure and yz projections of the "proper" (c) 
and (d) and complex non-collinear (CNC) (e) and (f) spiral 
or screw-type spin configurations are also shown. 



ccptcd spin-current model. This complex ground state 
is stabilized by the displacement of the oxygen atoms, 
which are also responsible for the multiferroic coupling 
predicted by Arima @. 

With a large S = 5/2 magnetic moment contributed 
by the Fe 3+ ions, CuFeO*2 (Fig. [Ha)) nas inspired great 
interest due to the magnetic frustration within each 
hexagonal plane [THJ] , which results in a four sublattice 
collinear spin configuration characterized by the prop- 
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FIG. 2: Energy versus momentum transfer taken from the 
triple-axis spectrometer at HFIR (a) and the time-of-fiight 
Cold Neutron Chopper Spectrometer (CNCS) at the SNS 
(b). (c) Predicted simulation of the INS data using intra- 
and interlayer interactions along with a lattice distortion and 
anisotropy. (d) Wavevector scans with different energy trans- 
fers for the CNCS data. 



agation wavevector [0.25,0.25,1-5]. A multiferroic state 
can be induced in this system by applying a magnetic 
field above about 7 T [lj| o r by doping with non- 
magnetic Al or Ga impurities [ill [l2| . The resulting 
multiferroic ground state is non-collinear and character- 
ized by the incommensurate propagation wavevectors Q 
= [H,H,1.5] where H 



0.20 and 0.30 \mUM- It is no- 
table that the spin excitation spectrum of the collinear 
non-multiferroic phase of pure CuFe02 softens at the in- 
commensurate wavevectors characteristic of the multifer- 
roic phase [l5[ . This softening can be simulated by low- 
ering the anisotropy energy [19j and the complete soften- 
ing of this mode signals the development of the multifer- 
roic non-collinear spin structure 

HUHlj]. The observed 
softening of the spin excitations of CuFe02 can be in- 
terpreted as a dynamical precursor of the multiferroic 
phase. 

Determining the spin configuration of multiferroic 
CuFe02 is nontrivial. A simple spiral structure, in which 
the average turn angle of Fe spins is 74°, would produce 
the observed peak at [0.20,0.20,1.5] but not the peak at 
[0.30,0.30,1.5] . Nakajima et al. [l]| proposed a "proper" 
spiral configuration (Fig. QJc) and (d)) with two alter- 
nating turn angles of 0° and 152° that will produce the 
two main observed peaks [H,H,1.5] (H 0.20 and 0.30). 

Recently, several groups have observed lattices distor- 
tions associated with displacements of the oxygen atoms 
(illustrated in Fig. [2(b)) EH, 23 1 . Incorporating these 
distortions, Fishman and Okamoto found that odd-order 
spin harmonics produce a magnetic ground state with a 
CNC spin configuration (Fig. [T|e) and (f)) [Hj]. The oxy- 



gen displacements produce a modulation in the nearest 
neighbor interactions J\ (Fig.QJb)) that breaks the equi- 
lateral symmetry (Fig. [2(b)) of the lattice with nearest- 
neighbor exchange interactions J± = j[ 2 ^ = J\ — K\/2 

(3) 

and J-J = J i +K±, where K\ is a measure of the atomic 
distortion. The distortion creates turn angles that fluc- 
tuate around 19° and 134°. For a distorted antifcrro- 
magnet, a CNC screw-type configuration is en erg etically 
favored over both simple and "proper" spirals [lj] . Elas- 
tic neutron scattering measurements cannot distinguish 
between the "proper" spiral and CNC configurations. 

The hexagonal lattice symmetry of CuFe02 provides a 
complex network of multiple intra- and interlayer super- 
exchange pathways (Fig. DIb)) [24| . The Hcisenberg 
Hamiltonian can be written as 



H 
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(i) 



where Si is the local moment on site i, D is the anisotropy 
energy, and the exchange coupling Jy between sites i and 
j is antiferromagnetic when Jy < 0. 

To incorporate the spin harmonics, we modified the 
classical approach described in Ref. [l4| by defining S z 
within any hexagonal plane as 

S Z {TL) = A{ cos(Q -x) + J2 c 2i+i cos[Q(2Z + 1) • x] 

i=i 

-J2 B ^i+i sin[(27r - Q)(2l + 1) • x}}, 

1=0 

(2) 

where the C2/+1 harmonics are produced by the 
anisotropy energy D and the B21+1 harmonics are pro- 
duced by the lattice distortion K\ . The square of these 
harmonics are proportional to the observed elastic inten- 
sities at odd multiples of Q and 2-7T — Q. The function 
S Z (R) is normalized so that the maximum of ^(R)! is 
S = 5/2. The perpendicular spin S y is given by 



where 



S„(R) = VS-5,(R) 2 8gn(g(R)), 



g(R) = sin(Q • x) + G x cos[(2tt - Q) 



(3) 



(4) 



and G\ is an additional variational parameter. 

Hexagonal spin planes were then stacked along the 
c axis and coupled by the exchange interactions J zn . 
The three-dimensional spin configuration was obtained 
by minimizing the energy on a large unit cell with length 
10 4 . The spin fluctuations about this three-dimensional 
spin configuration were evaluated using the method de- 
scribed in Ref. [25[ , where the frequencies and intensities 
of the SW excitations are evaluated simultaneously. We 
emphasize that a stable ground state is required to evalu- 
ate the spin dynamics: a 1/S spin- wave expansion cannot 
be performed starting from the "proper" helix sketched 
in Fig. DI.c). 
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FIG. 3: (a) Energy versus momentum transfer taken from 
the CNCS clearly showing the SW modes along the [0,0, L] 
direction with H = 0.2. The dispersion with L indicates the 
importance of interlayer interactions, (b) Predicted INS spec- 
tra using intra- and interlayer interactions along with a lattice 
distortion, (c) Intensity versus momentum transfer for vari- 
ous energies, (d) A three dimensional view of the Fe 3+ lattice 
structure illustrating the interlayer interactions. 



Single crystals of CuFco.965Gao.o3502 were grown using 
the floating zone technique. The crystals were oriented 
with two orthogonal wavevectors of [1,1,0] and [0,0,3] 
aligning in the horizontal plane. The INS measurements 
were carried out using the Cold Neutron Chopper Spec- 
trometer (CNCS) spectrometer at the Spallation Neutron 
Source (SNS) and the HB-1 triple-axis spectrometers at 
the High Flux Isotope Reactor (HFIR) at the Oak Ridge 
National Laboratory. The low energy excitations were 
measured at 1.4 K with final energy chosen at 5.1 meV 
at HB-1 and with incident neutron energy chosen at 3 
meV to ensure high resolution at CNCS. 

Figures HJa) and (b) present the INS data for the mul- 
tiferroic phase of 3.5% Ga doped CuFe02 along [H,H,0] 
with L = 1.5. "Gapless" excitations are observed at H w 
0.20 and 0.30 together with a "shoulder" at H w 0.08 and 
an intensity "hole" around H « 0.30 and E ps 1.0 meV. 
Figure [UJc) displays the the predicted INS spectra [26j |. 
which contains contributions from both the normal spin 
configuration with wavevector along [H,H,1.5] and the 
two twins with wavevectors along [H,0,1.5] and [0,H,1.5]. 
Strongly affected by the lattice distortion, the twin states 
account for the spectral weight around E m 1.5 meV and 
H « 0.30. The model also accurately produces the inten- 
sity "hole" at H w 0.30 as well as the "shoulder" at H « 
0.08, which is caused by the interlayer interactions. With 
the number of intricate features reproduced, we are con- 
fident that the dynamical "fingerprint" of doped CuFc02 
has been matched. 
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FIG. 4: (a) Intensity versus energy for the [0.20,0.20,1.5] 
(blue squares) and [0.30,0.30,1.5] (red circles) modes along 
the [H,H,0] direction with L = 1.5. The ratio of mode intensi- 
ties shows that the main [0.20,0.20,1.5] mode is three times as 
intense as the [0.30,0.30,1.5] mode. The simulated intensity is 
consistent with this intensity ratio, (b) SW velocities for the 
two modes along the [H,H,0] direction with L = 1.5 together 
with the predicted values (gray circles). 



To compare the experimental and theoretical results 
more closely, we examined the SW velocities dw/dq at 
low energy (Fig. H^b)). The data points are determined 
by the energy cuts plotted in Fig. [5{d). Since the SW 
velocities are functions of the exchange interactions, this 
comparison assures us that the correct exchange param- 
eters have been used. The exchange parameters J n and 
J zn are quite close [26] to those predicted in Ref. [24[ 
from fits to the SW spectra of the tT44- phase. The exper- 
imental results show no evidence of an energy gap present 
in the multifcrroic phase within the experimental energy 
resolution. This is consistent with our model, in which 
rotations of the spin plane about the c axis cost no energy 
and a SW gap is absent. However, the magnetostrictive 
energy that produces the observed confinement [l3| of 
the spins to the [1,1,0] plane may induce a small gap in 
the excitation spectrum. 

While four harmonics are predicted, only the B\ com- 
ponent produces sufficient scattering intensity compar- 
ing to the main mode at [0.20,0.20,1.5] to be seen by our 
INS measurements [26[ . While the main SW mode is as- 
sociated with the ordering wavevector Q, the secondary 
mode centering at [0.30,0.30,1.5] is produced by the B\ 
harmonic at 2-7T — Q. The overall magnitude of these 
modes are in agreement with our calculated values. For 
lower Ga doping or larger anisotropy, the higher harmon- 
ics will become more important and the distribution of 
turn angles will become more uniformly distributed. If 
D is increased from 0.01 to 0.04 meV (the critical value 
above which the tT44- phase is stable), then the ampli- 
tudes C3 and B3 rise by factors of 2 and 4, respectively. 

Similar to the collinear excitation observed [l5| in pure 
CuFe02, the dispersive excitation along the [0,0, L] di- 
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rcction through the magnetic Bragg peak at H sa 0.20 
contains a minimum at L = 1.5 (Figs. [3{a) and (b)). 
This underscores the importance of the interlayer inter- 
actions J zn (Fig. EJd)) and the three dimensional char- 
acter of this system. In comparison to pure CuFe02 , the 
magnetic interactions J Z 2 and J Z 3 are somewhat weak- 
ened [2^1 , which may be attributed to the disorder caused 
by Ga doping. Figure Etc) shows the experimental cuts 
in energy demonstrating the splitting of the Goldstone 
mode at L m 1.5. 

Figure a) compares the integrated scattering inten- 
sity originating from the main mode at H = 0.20 and 
the secondary mode at H = 0.30. The main mode is 
about three times as intense as the secondary one. The 
simulation with B\ = -0.52 yields an intensity ratio of 
3.7, in good agreement with the measured value. The 
only well-defined modes in our model are the Goldstone 
modes emerging from [0.20,0.20,1.5] and [0.30,0.30,1.5]. 
A continuum of excitations borders the Goldstone modes 
due to the complex and incommensurate nature of the 
magnetic structure. This contrasts with the relatively 
large gapped excitation and associated magnon-phonon 
hybridized electro-magnon observed in TbMnOa 
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According to Ref. the spin-driven multiferroic be- 
havior of CuFc02 arises from the uniform charge trans- 
fer through the metal-ligand hybridization in the pres- 
ence of spin-orbit coupling. Using a combination of INS 
measurements and theoretical modeling, we have demon- 
strated that lattice distortions play an important role 
in determining the magnetic ground state of Ga-doped 
CuFe02- With a slight modification of the exchange 
parameters from pure CuFe02 and by adjusting the 
anisotropy and distortion energies, we have achieved re- 
markable agreement between the observed and predicted 
dynamical "fingerprint" of the CNC multiferroic phase. 
This complex ground state provides an alternative way 
to realize multiferroic coupling, where displacements of 
the oxygen atoms severely distort the spin configuration 
and produce the electric polarization. Consequently, it 
is clear that the spin-current model 0-S] is not the only 
mechanism responsible for multiferroic behavior and that 
many other frustrated magnets with rhombohedral or 
hexagonal symmetries may also exhibit the same form 
of multiferroic coupling as doped CuFe02 [IH, [29| . The 
symmetry of these materials makes them candidates for 
exotic magnetoelectric control like the recently reported 
magnetic digital flop of ferroelectric domains [30( . 
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